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We analyze the decay B° K+ K~ K° using 383 million BB events collected by the BABAR 
detector at SLAC to extract CP violation parameter values over the Dalitz plot. Combining all 
K+K-K" events, we find Acp = -0.015 ± 0.077 ± 0.053 and P^ff = 0.352 ± 0.076 ± 0.026 rad, 
corresponding to a CP violation significance of 4.8a. A second solution near ■n/2 — Peg is disfavored 
with a significance of 4.5ct. We also report Acp and Pejj separately for decays to 0(lO2O)if°, 
/o(980)_ft:°, and K+R-R^' with m^+j^- > l.lGeV/c^. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



In the Standard Model (SM), the phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix [H is the sole source of CP violation in the quark sec- 
tor. Due to interference between decays with and without 
mixing, this phase yields observable time-dependent CP 
asymmetries in meson decays. In particular, signifi- 
cant CP asymmetries in 6 — > sss decays, such as B'^ 
K^K~K^ Q, are expected 0, Deviations from the 
predicted CP asymmetry behavior for B^ K^K^ 
are expected to depend weakly on Dalitz plot (DP) po- 
sition 0| . Since the h — > sss amplitude is dominated 
by loop contributions, heavy virtualparticles beyond the 
SM might contribute significantly [y, |7[ . This sensitivity 
motivates measurements of CP asymmetries in multiple 
b sss decays @, i, @, ^ . 

Previous measurements of CP asymmetries in B^ — > 
K^K^K^ have been performed separately for events 
with K^K^ invariant mass (tok+k-) in the (/) mass [llj 
region, and for events excluding the (j) region, neglecting 
interference effects among intermediate states [3, 8, 10]. 
In this Letter we describe a time-dependent DP analy- 
sis of K^K~K^ decay from which we extract the 
values of the CP violation parameters Acp and (3eff by 
taking into account the complex amplitudes describing 
the entire i?° and B° Dahtz plots. We first extract the 
values of the parameters of the amplitude model, and 
measure the average CP asymmetry in B^ K^K^ 
decay over the entire DP. Using this model, wc then 
measure the CP asymmetries for the (f)K^ and fviK^ de- 
cay channels, from a "low-mass" analysis of events with 
iTT-K+K- < l.lGeV/c^. Finally, we perform a "high- 
mass" analysis to determine the average CP asymmetry 
for events with rrix+K- > 1-1 GeV/c^. 

The data sample for this analysis was collected with 
the BaBAR detector [l3| at the PEP-II asymmetric-energy 
e+e" collider at SLAC. Approximately 383 x 10^ BB 
pairs recorded at the T(4S') resonance were used. 

We reconstruct K^K^K^ decays by combin- 

ing two oppositely-charged kaon candidates with a K'^ 
reconstructed as -> tt+tt" (BO^.,) [II, if" ^ Tr^Tr" 
(B^o,), or Kl (5°^)). Each tt^tt^ candidate is 

formed from two tt" — > 77 candidates. Each photon 
has > 50 MeV and transverse shower shape consistent 
with an electromagnetic shower. Both tt^ candidates sat- 
isfy 100 < TO-y^ < 155 MeV/c^ and yield an invariant mass 
rriT^OT^a in the range —20 < m^a^a — rrij^o < 30MeV/c^. 
A Kl candidate is defined by an unassociated energy de- 



posit in the electromagnetic calorimeter or an isolated 
signal in the Instrumented Flux Return 

For each fully reconstructed meson (Bcp), we use 
the remaining tracks in the event to reconstruct the decay 
vertex of the other B meson (i?tag), and to identify its 
flavor (/tag [4] . For each event we calculate the difference 
At = tcp — ttag between the proper decay times of the 
Bcp and Stag mesons, and its uncertainty a^t- 

We characterize S°+_) and -8°^^ candidates using two 

kinematic variables: the beam-energy-substituted mass 

Toes and the energy difference AE [8] . The signal region 

(SR) is defined as toes > 5.26 GeV/c^, and |A£:| < 0.06 

GeV for or -0.120 < A£: < 0.06 GeV for 

For B°, the SR is defined by -0.01 < AE < 
I— I i^) ^ 

0.03 GeV |8|, and the missing momentum for the entire 
event is required to be consistent with the calculated 
laboratory momentum. 

The main source of background is continuum e^e^ — > 
qq {q = u, d, s, c) events. We use event-shape variables to 
exploit the jet-like structure of these events in order to 
remove much of this background . 

We perform an unbinned maximum likelihood fit to 
the selected K^K^K^ events using the likelihood func- 
tion defined in Ref. [s*] . The probability density function 
(PDF), Vi, is given by 

Vr = V{mEs) ■ V{AE) ■ Vlow (1) 

• 'PDp{mK+K-,COSeH,At,qta.g) (g) TZ{At , (T At) , 

where i = (signal, continuum, BB background), and TZ 
is the At resolution function [4]. For B'^^^, V{m^s) is not 
used. Vhovj is a PDF used only in the low-mass fit, which 
depends on the event-shape variables and, for B'^^^ only, 
the missing momentum in the event Q . We characterize 
5° [B^) events on the DP in terms of rrix+x- and cos 9h, 
the cosine of the helicity angle between the {K~) and 
the K° (F°) in the rest frame of the K+K' system. The 
DP PDF for signal events is 

VDP^dT-eimK+K-,cos9H)-\Jl (2) 

where dT is the time- and fiavor-dependent decay rate 
over the DP, e is the eflaciency, and J is the Jacobian of 
the transformation to our choice of DP coordinates. 
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The time- and flavor-dependent decay rate is 

At|/r 



dr e 



i^r + i^i (3) 

+ gtag 2Im {^AA*) sin AnidAt 



qta 



A\ J cos AnidAt 



where r and Auid are the Ufetime and mixing frequency 
of the -6° meson, respectively The parameter ^ — 

??cpe"^''^, where /3 = arff(-Vc<iV;y FtdV^;) and Vg,' are 
CKM matrix elements The CP eigenvalue tjcp — 
1 (-1) for the -RTg mode. We define the amphtude 



A for decay as a sum of isobar amplitudes [14| . 



A {rriK+K- , cos6h) — ^ ^ r 

r 

= ^c,(l T6r)e^^'^'-^''-^ • /.(mK+if-,cos0H) 



(4) 



where the minus signs are associated with the A, the 
parameters and tpr are the magnitude and phase of 
the amplitude of component r, and we allow for differ- 
ent isobar coefficients for i?° and decays through the 
asymmetry parameters br and Sr- 

Our isobar model includes resonant amplitudes 0, /o, 
Xco(l-P)j and Xq{1550) [iMIlBl; non-resonant terms; and 
incoherent terms for decay to D^K^ and DjK^. 
For each resonant term, the function fr — x Tr x Zr 
describes the dynamical properties, where Fr is the Blatt- 
Weisskopf centrifugal barrier factor for the resonance de- 
cay vertex 17[, is the resonant mass-lineshape, and 
Zr describes the angular distribution in the decay 
The barrier factor F,. = 1/^1 + [17] for the 0, 

where q is the momentum in the 4> rest frame and 
R = 1.5 GeV~^; Fr = I for the scalar resonances. For 
(j> decay Zr ^ q ■ p, where p is the momentum of the 
in the (j) rest frame, while Z,. = 1 for the scalar decays. 
We describe the (j), Xn(1550), and Xcoi^P) with relativis- 
tic Breit-Wigner hneshapes [l4|. For the (j) and Xco(l^) 
parameters we use average measurements [l^ . For the 
Xo(1550) resonance, we use parameters from our analysis 
of the -6+ K^K^K^ decay [15]. The /o resonance is 
described by a coupled-channel amplitude [l^, with the 
parameter values of Ref . [l^] . 

We include three non-resonant (NR) amplitudes pa- 
rameterized as fNR,k = exp(— am^), where the param- 
eter a = 0.14 ± 0.01 c^/GeV'^ is taken from measure- 
ments of 5+ — > K^K^K^ decays with larger signal 
samples [isl. We include a complex isobar coefficient 
for each component k = {K+ R- , K+ , R- K°). 

PDFs for qq background in i?" R^R^ R2. are mod- 
eled using events in the region 5.2 < mEs < 5.26GeV/c^. 
The region 0.02 < A£; < 0.04 GeV is used for 5°^,. Sim- 
ulated BB events are used to define BB background 
PDFs. We use two-dimensional histogram PDFs to 
model the DP distributions for qq and BB backgrounds. 



We compute the CP asymmetry parameters for com- 
ponent r from the asymmetries in amplitude {br) and 
phase (6r) given in Eq. The rate asymmetry is 



Acp,r 



I Ar \ \ Ar 
\A.r I ] 



-2br 



(5) 



and Peff,r = (3 + Sr is the phase asymmetry. 

The selection criteria yield 3266 1611 Bj'oo,, and 

27513 candidates which we fit to obtain the event 
yields, the isobar coefficients of the DP model, and the 
CP asymmetry parameters averaged over the DP. The 
parameters br and Sr are constrained to be the same for 
all model components, so in this case Acp.r — Acp and 
Peff^r = Peff- Wc find 947 ± 37 144 ± 17 fiO^,, 

and 770 ± 71 B^^^ signal events. Isobar coefficients and 
fractions are reported in Table [H and CP asymmetry re- 
sults are summarized in Table HH The fraction Tr for 
resonance r is computed as in Ref. [l^. Note that there 
is a ±7rrad ambiguity in the Xco{^P)K^ phase. 



TABLE I: The isobar amplitudes Cr, phases </?r, and fractions 
Tr from the fit to the full K+ K' K" DP. The three NR com- 
ponents are combined for the fraction calculation. Errors are 
statistical only. Because of interference, X^.?> 7^ 100%. 



Isobar Mode 


Amplitude Cr Phase ipr ( rad) J?> (%) 


Xo(1550)ii'° 
{K+K-)nrK° 
iK+K°)NRK- 
{K-K°)nrK+ 


0.0085 ± 0.0010 -0.016 ±0.234 12.5 ±1.3 
0.622 ±0.046 -0.14 ±0.14 40.2 ± 9.6 
0.114 ±0.018 -0.47 ±0.20 4.1 ± 1.3 
1 (fixed) (fixed) 
0.33 ± 0.07 1.95 ± 0.27 112.0 ± 14.9 
0.31 ±0.08 -1.34 ±0.37 


D-K+ 
D7K+ 


0.0306 ± 0.0049 .2:33 ± 0.54 3.0 ±1.2 
1.11 ±0.17 3.6 ±1.5 
0.76 ±0.14 1.8 ±0.6 



In Fig. [U we plot twice the change in the negative 
logarithm of the likelihood as a function of f5eff ■ We find 
that the CP-conserving case of fieff = is excluded at 
4.8(7 (5.1(t), including statistical and systematic errors 
(statistical errors only). Also, the interference between 
CP-even and CP-odd amplitudes leads to the exclusion 
of the Pe^ff solution near n/2 — /? at 4.5cr (4.6cr). 




Statistical en'ors only 
Statistical and systematic errors 



0.6 0.8 



1.2 



1.4 1.6 



FIG. I: The change in twice the negative log likelihood as a 
function of jS^ff for the fit to the whole DP. 
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FIG. 2: The distributions of niK+K- for signal- weighted [23| 
data in (a) the entire DP and (b) the low-mass region. 
Insets show distributions of cos 9h ■ The histograms are pro- 
jections of the fit function for the corresponding result. 
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1 ] 



-8 -6 -4 -2 2 4 6 8 

At (ps) 

FIG. 3: The raw asymmetry between S"- and B^-tagged 
signal- weighted [23| events for B'^^ in (a) the low-mass re- 
gion and (b) the high-mass region. The curves axe projections 
of the corresponding fit results. 



We also measure CP asymmetry parameters for events 
with TTiK+K- < 1-1 GeV/c^. In this region, we find 1359 
348 and 7481 5°^, candidates. The fit yields 

282 ± 20, 37 ± 9 and 266 ± 36 signal events, respectively. 
The most significant contributions in this region are from 
(j>K'^ and f^K^ decays, with a smaller contribution from 
the low-mass tail from non-resonant decays. In this fit 
we vary the amplitude asymmetries hr and 5r for the (j) 
and /o, while the other components are fixed to the SM 
expectations of fieff = 0.370 rad and Acp — [2l|. We 
also vary the isobar coefficient for the 0, while fixing the 
others to the results from the whole DP fit. There are two 
solutions with likelihood difference of only A log L — 0.1. 
Solution (1) is consistent with the SM, while in Solu- 
tion (2) (ieff for the /o differs significantly from the SM 
value (Table [11} • The solutions also differ significantly 
in the values of the isobar coefficient. There is also a 
mathematical ambiguity of ±7r rad on (3eff for the (f>, with 



a corresponding change of ±7rrad in the solution for ip^. 
This ambiguity is present for both solutions. The fit cor- 



22| 



and 20032 
withmx+K- > 1-1 GeV/c^ 



2384 B°^^^ 



1406 



relation between the (f> and /o in 5r is 0.71 

Finally, we perform a fit to extract the average CP 
asymmetry parameters in the high-mass region. In the 

selected events 
we find signal yields of 673 ± 
31, 87 ± 14 and 462 ± 56 events, respectively; the CP 
asymmetry results are shown in Table [III We find that 
for this fit the CP-conserving case of (3eff = is excluded 
at 5.1a, including statistical and systematic errors. 

Figure [2] shows distributions of the DP variables 
rriK+K- and cos 9h obtained using the method described 
in [23[. Figure [3] shows the At-dependent asymmetry be- 
tween B^- and S^-tagged events. 



TABLE II: The CP-asymmetries for B° K+ K' K° for the 
entire DP, in the high-mass region, and for i^Tf " and foK^ in 
the low-mass region. The first errors are statistical and the 
second are systematic. The solutions (1) and (2) from the 
low-mass fit are discussed in the text. 





Acp 


f3eff ( rad) 


Whole DP 


-0.015 ±0.077 ±0.053 


0.352 ± 0.076 ± 0.026 


High-mass 


-0.054 ±0.102 ±0.060 


0.436 ± 0.087 1 o.oM 


(1) <t>K" 
(1) /oi^° 


-0.08 ±0.18 ±0.04 
0.41 ±0.23 ±0.07 


0.11 ±0.14 ±0.06 
0.14 ±0.15 ±0.05 


(2) ^K- 
(2) /oT^o 


-0.11 ±0.18 
-0.20 ±0.31 


0.10 ±0.13 
3.09 ±0.19 



Systematic errors on the CP-asymmetry parameters 
are listed in Table IIIII The fit bias uncertainty includes 
effects of detector resolution and possible correlations 
among the fit variables determined from full-detector 
simulations. We also account for uncertainties due to 
the isobar model: experimental precision of resonance pa- 
rameter values; alternate Xo( 1550) parameter values [16l |: 
and, in the low- and high-mass fits, the statistical uncer- 
tainties on the isobar coefficients determined in the fit 
to the whole DP. Other uncertainties common to many 
BaBAR time-dependent analyses, including those due to 
fixed PDF parameters, and possible CP asymmetries in 
the BB background are also taken into account [1, [ij] ■ 
Uncertainties due to fixed PDF parameters are evaluated 
by shifting the fixed parameters and refitting the data. 
As a cross-check, we perform the analysis using B^^_^ 
alone and find results consistent with those in Table |TT1 



TABLE HI: A summary of the systematic errors on the CP 
asymmetry parameter values. 



Source 


Whole DP 

Acp Peff 


High-mass 
Acp Peff 


(l>K^ foK" 

Acp Peff Acp Pejf 


Fit Bias 
Isobar model 
Other 


0.003 0.001 
0.004 0.009 
0.052 0.024 


0.014 0.008 
025 

U.UZO _o.024 

0.053 0.018 


0.03 0.06 0.06 0.03 
0.00 0.01 0.01 0.03 
0.02 0.01 0.03 0.02 


Total 


0.053 0.026 


060 

U.UUU _o.031 


0.04 0.06 0.07 0.05 
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In summary, in a sample of 383 x 10^ BB meson 
pairs we simultaneously analyze the DP distribution and 
measure the time-dependent CP asymmetries for B'^ 
K'^K~K^ decays. The values of jSeff and Acp are consis- 
tent with the SM expectations of /9 ~ 0.370 rad, Acp — 
[2l|. The signficance of CP violation is 4.8tT, and we 
reject the solution near 7r/2 — /? at 4.5cr. We also mea- 
sure CP asymmetries for the decays B^ — > (f)K^ and 
B^ JqK'-^ where wc find f3eff lower than the SM expec- 
tation by about 2a. The CP parameters in the high-mass 
region are compatible with SM expectations, and we ob- 
serve CP violation at the level of 5.1a. 
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